Drifting instabilities of cavity solitons in vertical cavity surface-emitting lasers 
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In this paper we study the formation and dynamics of self-propelled cavity solitons (CSs) in a 
model for vertical cavity surface-emitting lasers (VCSELs) subjected to external frequency selective 
feedback (FSF), and build their bifurcation diagram for the case where carrier dynamics is elimi- 
nated. For low pump currents, we find that they emerge from the modulational instability point of 
the trivial solution, where traveling waves with a critical wavenumber are formed. For large currents, 
the branch of self-propelled solitons merges with the branch of resting solitons via a pitchfork bifur- 
cation. We also show that a feedback phase variation of ±27r can transform a CS (whether resting 
or moving) into a different one associated to an adjacent longitudinal external cavity mode. Finally, 
we investigate the influence of the carrier dynamics, relevant for VCSELs. We find and analyze 
qualitative changes in the stability properties of resting CSs when increasing the carrier relaxation 
time. In addition to a drifting instability of resting CSs, a new kind of instability appears for certain 
ranges of carrier lifetime, leading to a swinging motion of the CS center position. Furthermore, for 
carrier relaxation times typical of VCSELs the system can display multistability of CSs. 
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I. INTRODUCTION 

In recent years jthere has been significant progress in 
the experimental [ll, and theoretical 0,13 study of self- 
localized states in vertical cavity surface-emitting lasers 
(VCSELs) subject to frequency selective external optical 
feedback (FSF). These systems are attractive from an 
experimental point of view because they can be imple- 
mented with basically off-the-shelf optical components, 
and they do not require an optical holding beam to sup- 
port self-localized solutions, also known in this context as 
cavity solitons (CSs). These systems are known as Cav- 
ity Soliton Lasers, because for certain parameter regions 
laser emission takes place only in localized structures. 
An important property of this system is its invariance 
under phase transformations. Hence each CS is a mem- 
ber of a continuous family of phase-equivalent localized 
structures. Further, for the same parameters there may 
be several families of CSs with different frequencies. In 
contrast, localized structures in driven-cavity systems are 
phase- and frequency-locked to the holding beam. 

There are many possible applications of CSs but 
we will emphasize here only the most relevant one for 
the contents of this article, namely an optical delay line 
or shift register. This element is essential in communi- 
cation systems for delaying pulses when routing several 
sequences of bits. This function was proposed to be im- 
plemented on the basis of spatially drifting CSs [1] , where 
the drift of CSs was caused by a gradient of the system 



parameters. Thus, one can inject a sequence of pulses 
at one transverse position, and read out a copy of this 
sequence at another position at some later time (delay 
time of the shift register) . 

In contrast with gradient-induced CS motion, the pos- 
sibility of self-motion was shown in a holding beam sys- 
tem [7| with thermal effects. Spontaneous motion of vor- 
tices in lasers and laser amplifiers with saturable absor 
tion has been also studied in [8|. In a previous work 
we have shown the existence of 2D self-propelled CSs in 
a VCSEL with FSF, even without thermal effects. This 
is an alternative mechanism to induce motion in delay 
lines. Since no parameter gradient has to be imposed 
and controlled, there is, in principle, no limit to the drift 
distance, and hence delay time of the shift register. Since 
CSs moving in the transverse direction can be reflected 
by the boundaries of the device 0] the delay line can be 
folded, enabling longer paths for a given VCSEL aper- 
ture, and therefore improving the delay-bandwidth prod- 
uct. 

In this article we analyze in detail the bifurcation di- 
agram of self-moving solitons elucidating the conditions 
for the appearance of motion in the ID case. In com- 
parison with [3] here we address the role of the feedback 
phase, which plays an important role in determining the 
longitudinal- mode frequency of the CS. For a given feed- 
back phase several CS solutions with different frequen- 
cies (associated to adjacent external cavity modes) are 
found to coexist. For each solution, sweeping the feed- 
back phase allows to change continously the frequency 
of the CS. Applying a 2tt sweep in the feedback phase 
shows a smooth transition from one CS solution branch 
to the adjacent one. The number of coexisting solutions 
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increases with the delay time. We also investigate the 
role of the carrier dynamics, which was neglected in [J, 
and analyze its influence on the stability of CSs. In par- 
ticular we show the existence of a swinging instability in 
which the position of the CS maximum starts to perform 
growing oscillations near the initial position. 

The paper is organized as follows: In Section|ll]we dis- 
cuss the system and the models considered. In Section HIIl 
we describe the properties of self-propelled CSs and for- 
mulate the equations for their semianalytical calculation 
using a Newton method. In this section we present also 
the bifurcation diagrams of moving and resting solitons 
explaining their connection with the stability properties 
of the nonlasing background solution and discussing the 
relation between the ID and the 2D problems. In section 
IIVI we study the possibilities of controlling the solitonic 
states by varying the feedback phase and show the co- 
existence of multiple states, both moving and resting. In 
Section |V] we study the influence of the carrier dynamics, 
leading to the swinging instability and to the stabiliza- 
tion of resting solitons. Finally, in Section IVII we give 
some concluding remarks. 



II. SYSTEM AND MODEL 

Following [l|, 0, [1, 3 , we consider the set-up sketched 
in Fig. [T]), where light is fed back to the VCSEL after 
being filtered in frequency. The polarizer allows feedback 
only for the linearly polarized field component that is 
dominant in the solitary laser, so that the field can be 
treated as a scalar. 
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FIG. 1: (Color online). Scheme of the system. Fi and F2 are 
the focal lengths of the first and second lenses respectively. 



The dynamics of the complex envelope of the electrical 
field E{x,t), the filtered feedback field F{x,t) and the 
carrier number N(x, t) can be described by 



' 1^ = k(1 + ia){N - 1)E - iAj_E + F + iuj.E, 



-ae"t'XE{t-T), 



(1) 



where k is the decay rate of the field in the cavity; a 
is the linewidth enhancement factor describing phase- 
amplitude coupling; fi is the pump current normalized 
to be 1 at the threshold of the solitary laser; Aj_ = ^3 
describes diffraction and 7 = ^ is the carrier relaxation 
rate. The filter is considered to have a Lorentzian fre- 
quency response p^ . [Til . [T^ . and its central frequency 
is taken as referenceHence lOs the frequency of the axial 
mode of the solitary laser at threshold {N = 1). The 
other filter parameters are cr, the feedback strength, cf) 
the feedback phase, r the delay time in the feedback loop, 
and A the filter bandwidth. 

Eliminating the carriers adiabatically model ([1]) re- 
duces to 



dE 
dt 



dF 
dt 



+ ia)E + k(1 + ia)Yf^ 
-iAj_E + F + iujsE, 

-\F + ae"f'XE{t-T), 



(2) 



For a fixed feedback phase 0=0 this was the model 
considered in Q. We will use this simplified model in 
Sect. lIIII to determine the bifurcation diagrams of moving 
and resting solitons as well as in Sect. IIVI to discuss the 
role of the feedback phase. While the simplified model 
allows for easier calculation, the typical relaxation times 
of the system are in fact « 0.01 ns, « 0.2 ns 
and Ti « Ins (throughout this paper time is measured 
in nanoseconds). Since the carrier density is the slowest 
variable, adiabatic elimination of the carriers is not fully 
justified. Nevertheless, as discussed in Sect. |V] where 
we address the infiuence of the carrier dynamics by con- 
sidering the full model ((1]), the reduced model correctly 
predicts the existence of moving and resting stable CSs, 
albeit with slight changes in the parameter regions where 
they are observed. 



III. RESTING AND MOVING CAVITY 
SOLITONS 

We consider in this section the simplified model ^ 
for feedback phase (f) — 0, as in Q. Fig. [2] shows the 
marginal stability curves for the trivial solution E ^ 0. 
This solution is stable for small pump values (below line 
D in the figure). At ^ = /is = 1 — ct/k = 0.4 it becomes 
unstable to perturbations with the frequency of the fil- 
ter maximum and with a finite transverse wavenumber. 
For values of the pump between lines D and B the for- 
mation of a complex spatiotemporal regime is observed. 
For pump values in between line B and the thick line at 
/i w 1 the trivial solution is stable again. In this region, 
which we will call gap-region, we observe the spontaneous 
formation of self-localized states, both resting CSs and 
self-propelled CSs [|. In the ID case, we observe the 
formation of self-propelled CSs for almost all the values 
of the pump in the gap-region whereas in the 2D case this 
region splits in 3 subregions: a) close to B, where moving 
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FIG. 2: (Color online). Marginal stability curves of the non- 
lasing state in the VCSEL with external feedback. The inset 
shows the frequencies of the modes that become unstable be- 
tween lines D and B in the main figure. The gray (cyan) 
curve corresponds to the boundary of excitation of external 
cavity modes with frequencies close to tj = (the filter maxi- 
mum) . The line A corresponds to value of the pump at which 
the marginal instability curve gray (cyan) intersects the axis 
fci = 0. The thin dark (red) curve corresponds to the bound- 
ary of excitation of the two modes at either side of the domi- 
nant mode with frequencies lo ~ ~ ±1.25 (the boundary 
is practically the same for both modes for the parameters con- 
sidered here). The thick black line corresponds to the excita- 
tion of modes with frequencies close to those of the solitary 
laser. The circles represent the square of the wavenumber 
of the far-field maximum for self-propelled ID CSs. Stable 
localized solutions can exist in the region between the line B 
and the thick black line, in which the trivial solution is stable. 
Parameters: = 0, cr = 60, q = 5, A = 2.71, t = 5, u>s = 250, 
K = 100. 



2D CSs arc excited, b) intermediate values, where resting 
2D CSs are observed, and, c) values closer to 1, where no 
solitons can be excited 

The system is invariant under translational and global 
phase transformations, and thus the operations E(x) — > 
E{x + xq) and E{x) E{x)e^'^'^' transform solutions 
in solutions. So there is a whole family of CSs which 
have the same distribution of the far field absolute value 
\E{k)\. The members of the family are identified by two 
continuous parameters, their location and global phase. 
Moving solutions break the left-right symmetry of the 
system so that two equivalent families of self-propelled 
CSs exist, one with a positive k^ax moving to the left 
and its specular image with negative k^ax and moving 
to the right. 

The spatial distribution of a self-propelled ID CS for 
^ = 0.65 is shown in Fig. O The near field is character- 
ized by an exponential spatial localization. The far field 
shows a typical cascade of energy from small to large 
wavenumbers due to nonlinearity. The far field is cen- 
tered off-axis (at kmax = —2.54 in Fig. [3]) and is asym- 



1 r 

0.8 
gO.6 
y-0.4 

0.2 



i 



40 r 
_30 

sio 

10 



0.2 0.4 0.6 0.8 

(c) 



-200 200 
(d) 




FIG. 3: Shape of moving ID CSs; (a) Near field, (b) far field, 
(c) near field in semilogarithmic scale, and (d) far field in 
semilogarithmic scale. Parameters are the same as in Fig. |2] 
and iJ. — 0.65. 



metric: we note that the part of the far field to the right 
of the maximum is smoother than the left part. Self- 
propelled CSs move at a constant velocity, which is pro- 
portional to —kjnax- 

The residual level at ~ 10^^^, both in the near and far 
fields, is indicative of the numerical precision. The spatial 
profile of the soliton shown in Fig. [3] has been obtained 
from numerical integration of Eq. ([2]) where the initial 
condition breaks the left-right degeneracy. As for the nu- 
merical integration, we should note that accurate results 
arc difiicult to obtain using usual finite difference meth- 
ods, because large k± transverse wave vectors lead to the 
presence of very high frequencies, and hence they require 
very small temporal steps. For the numerical integration 
we have used the two-step pseudo-spectral method de- 
scribed in the appendix of [l3|, where in Fourier space 
the linear part of the equations is handled analytically, 
thus avoiding the problem with the high frequencies asso- 
ciated to high wavenumbers due to diffraction. Nonlinear 
terms are evaluated in real space. The pseudo-spectral 
method can be easily adapted to equations with delay 
such as the ones considered here. We have taken the 
temporal stepsize much smaller than all the character- 
istic times of the system and we have verified that the 
selected stepsize allows to reproduce correctly the thresh- 
old of the background instability predicted by the linear 
stability analysis, and the threshold of the drifting insta- 
bility obtained by the Newton method described below. 

One of the main objectives of this work is, however, to 
compute the whole branch of self-propelled states, includ- 
ing the region where they are unstable, hence we need to 
implement a Newton method, similar to what was done 
for resting 2D solitons in for moving solutions. We 
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will seek solutions of the form: 

" E{x,t) = Eo{x + vt)e^'^\ 

F{x,t) = Fo{x + vt)e''^* , 
Substituting ^ into ^ we have the following equations 



(3) 



f -^v + i-nil + ia) 



iuj)Eo 



+ ia)^ 



1+1 Bo 



+ i 



Fo = 0, 



[ -f^v - iujFn - AFo + ae^'^XEo{x' - vT)e'^^^'' = 0, 

(4) 

where x' = x + vt\s the spatial coordinate in the moving 
reference frame. In Fourier space 



Fo{k')=T{Fo{x')), 
Eo{k') = TiEoix')), 



T[Eo{x' -vT))=^^'^^Eo{k'), 



(5) 
(6) 
(7) 



the set of differential equations (j4|) is transformed into a 
set of algebraic equations. The second equation is linear 
and can be solved directly 



Fo{k') 



ere 



i(p '^^i{k' v—Lij)i 



■Eo{k'), 



X + i{LU + k'v) 
leading to a single equation for the field Eq 

T-\E^{k')) 



AEo{k') + BT 



l + \T-^{Eo{k'W 



(8) 



0, (9) 



where 



A = 



B = 



— k(1 + ia) + i{uLis — w — k'v — fc'^) 



\ + i{bj + k'v) ' 
k{\ + ZQ;)/i. 



(10) 

(11) 

(12) 



Equation ^ is analogous to the Eq. (3) of [j| for station- 
ary solutions, while it also supports more complex solu- 
tions, like self-propelled CSs. In practice Eq. ([9]) is solved 
numerically by discrctizing k' in n values which leads to 
a set of 2n coupled equations. There are 2ri,4-2 unknowns 
[£'0(^0 (* = 1, .-.n), V and to] and two constraints asso- 
ciated to the translational and global phase invariance. 
This set can be solved using a Newton method start- 
ing from an initial guess prepared by direct integration 
of ([21). At each Newton iteration, imposing that in the 
moving frame the location of the CS does not change 
allows for the determination of v. Similarly imposing 
that the global phase does not change allows for the de- 
termination of Lj. Equivalent solutions can be obtained 
by applying a translation or a global phase change. The 
specular family of self-propelled CSs can be obtained just 
by reflection. Furthermore, once a precise solution for a 
given set of parameters is obtained, it is possible to build 
the full branch of ID self-propelled CSs for different pa- 
rameter values using continuation methods. 



The branch of ID self-propelled CSs is displayed in 
Fig. [5] as circles indicating the square of the wavenum- 
ber of the far-field maximum kmax- A bifurcation di- 
agram showing the maximum amplitude of resting and 
self-propelled CSs is shown in Fig. |4l The dependence 
of the frequency lo, square of the velocity and spatial 
width of ID self-propelled CSs on the pump current is 
shown in Fig. [51 On one side the branch starts from 




FIG. 4: (Color online). Bifurcation diagram of the homoge- 
neous solutions (dotted green curve), ID-resting CSs (dash- 
dotted red curve), 2D-resting CSs (solid blue curve), and ID 
self-propelled CSs (circles). |i5,niiii| corresponds to the maxi- 
mum field amplitude for localized structures and to the field 
amplitude for the homogeneous solution. The vertical dash- 
dotted (blue) lines are the boundaries of the stable and un- 
stable regions of the trivial solution corresponding to B and 
D and = 1 in Fig. [2l Parameters are the same as in Fig. [2] 

pump current ^ = fijj, where self-propelled CSs appear 
with zero amplitude (Fig.^l) and with a kmax which coin- 
cides with the critical transverse wavenumber of the zero 
solution (intersection of line D and gray (cyan) curve in 
Fig. 121). Furthermore, sclf-propcUcd CSs originate with 
a negative detuning with respect to the filter central fre- 
quency (See Fig. inija)), with a finite velocity (Fig. jS^b)) 
and with an infinite width (Fig. EJc)). This last char- 
acteristic is similar to the appearance of resting CSs in 
points A and B of Fig. 21 (see also Q])- The shape of self- 
propelled CSs close to point D is displayed in Fig. |6l for 
/i = 0.405. The real part of the field performs many oscil- 
lations within a HWHM of a sech-like envelope. The fact 
that the width diverges while the oscillation wavcvcctor 
remains finite (corresponding to the critical wavenum- 
ber of the zero solution) does not allow to follow the 
branch close to bifurcation point just by rescaling the 
spatial length scale. It requires increasing of number of 
discretization points (2048 points were used in our ID 
calculations). As might be expected, it can be proved 
analytically that in this limit the moving CS is an en- 
velope soliton of the critical traveling-wave solution, and 
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FIG. 5: (Color online). Dependencies of the frequency (a) , the 
square of the velocity (b) and the half width at half maximum 
(HWIfM) (c) of the ID self-propelled CS on the pump current. 
The vertical lines D and M' mark the starting point and the 
end of the branch respectively. 

moves with its group velocity. 

The moving-CS branch in Fig. [5]) defines only exis- 
tence, and says nothing about stability. It is clear that 
the entire section of the branch below line B must corre- 
spond to unstable CS, because of instability of the zero- 
amplitude background state. Above line B, we find nu- 
merically that the moving cavity solitons are usually sta- 
ble, though those associated with sidebands may show 
instability (sec section HV|) . 

(a) (b) 
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FIG. 6: Shape of a moving ID CS close to the bifurcation 
point D of Fig. 3] (a) Near field amplitude (solid thick line) 
and its real part (solid thin line); (b) far field amplitude, 
/i = 0.405. 

It is worthwhile to emphasize that the motion of local- 
ized states is connected with the modulational instability 
of the background leading to travelling waves with a fi- 
nite critical transverse wavevector. This suggests a way 
of eliminating the drifting instability of resting solitons 
by removing the finite wavelength instability. This can 
be done by varying the detuning w^, which shifts the crit- 
ical transverse wavevector along the fc^ axis (see H). If 



instead we want to enhance the self-propelled solitons (in- 
crease the velocity) we need, then, to increase the critical 
wavenumber. 

Increasing the pump, the value of |fc^„a3.| and the as- 
sociated value of decrease (see Fig. [5] and Fig. [5ljb)). 
Simultaneously, the self-propelled CS becomes narrower 
(Fig. m (c)), taller (Fig. and its frequency approaches 
the one of the filter (Fig. [5] (a)). 

The branch of self-propelled CSs ends at point Af, 
where the moving solution merges with the resting one 
(Fig. [4]). In fact, at A/' a drift pitchfork bifurcation takes 
place. Decreasing ^, two stable self-propelled CSs with 
velocities v and —v originate from the resting CS that 
becomes unstable. The squared variables, and femaKi 
scale linearly from zero with the distance to the bifur- 
cation point (see Figs. Ob) and [2]), as expected for a 
supercritical pitchfork. 

Finally, some remarks on the 2D case arc in order. In 
2D, self-propelled CSs can move in any arbitrary direc- 
tion. The drift is associated to an asymetric cross-section 
along the the direction of motion while the section in the 
orthogonal (transverse) direction is symmetric. Due to 
the lack of rotational symmetry, finding the branch of 
self-propelled 2D CSs using a Newton method as above 
requires solving a full 2D problem. This is too demanding 
computationally for us to obtain a bifurcation diagram 
to compare with the ID case. Some 2D results may be 
inferred, however. Beyond the resting and self-propelled 
ID CSs, figure Fig.[3]also shows the resting 2D CS branch 
(solid (blue) curve) and the homogeneous solutions (dot- 
ted (green) curves) from Comparing the curves for 
the ID and 2D cases one may expect that the 2D self- 
propelled CS branch will start at the same point Z), with 
zero amplitude, infinite width, finite velocity and wave- 
length equal to the critical transverse wavenumber of the 
zero solution. Increasing the pump, self-propelled 2D 
CSs will become narrower and their speed will decrease. 
The branch of self-propelled 2D CS ends at point M in 
Fig. m where it merges with the resting CS branch. The 
point M has been obtained from numerical simulations 
of eqs. ^ for the 2D case [3]. The natural extension of 
the drift pitchfork bifurcation observed in ID to 2D is a 
drift circle pitchfork bifurcation [l^ in which a CS mov- 
ing in an arbitrary direction originates while the resting 
CS becomes unstable. 



IV. CS CONTROL VIA FEEDBACK PHASE 

In the previous sections we have discussed the proper- 
ties of model ([2]) with a fixed feedback phase. However, 
in delayed systems, the variation of this parameter leads 
to transitions between adjacent external cavity modes 
0, ■ Therefore we consider here the influence of feed- 
back phase (/) on the solitonic branches in model ([2]). In 
Fig. [7] the solid lines show the branches of resting ID 
CSs as a function of 0. The branches have an elliptical 
shape. The part of the branch with larger amplitude, and 
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lower frequency, is associated to external cavity modes, 
while the o ppo site part is associated to anti-modes, sad- 
dle points [131 that act as separatrices in phase-space. 
Consider, for instance, the branch shown as a thick line. 
At = the higher amplitude CS Z\ is stable while 
the lower amplitude CS Z2 is a saddle. As the feedback 
phase is varied, the amplitude and the frequency of the 
CS changes continuously. The branch covers a feedback 
phase interval larger than 47r. Since the feedback phase 
is 27r periodic, for some values of the feedback phase six 
solutions exist, while for others there are only four, as is 
clear from Fig. (T^. 

(a) 
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FIG. 7: (Color online). Bifurcation diagram of ID CSs as a 
function of the feedback phase 0: (a) Maximum amplitude; 
(b) frequency. The solid thick curve (blue online) corresponds 
to resting CSs; the open (filled) circles correspond to the un- 
stable (stable) self-propelled CSs. The vertical dashed lines 
correspond to = 0, (/> = ±27r. From the 27r periodicity, 
translating the thick curves gives the set of equivalent gray 
curves. Here [i, = 0.6 and the other parameters as in Fig. (2] 

For example, for = 0, resting CSs Z3, Z4, Z5, Zg exist 
in addition to Z\ and Zi. The coexisting "mode" states 
(Zi, Z3, Z5) have significantly different frequencies, with 
a separation approximately of 27r/T w 1.25, similar to the 
frequency separation of the marginal instability curves in 
Fig. [21 The coexisting antimodes for a given phase are 
also associated to neighbor external cavity antimodes (Z2 
to the central, Z4 to the high frequency and Zg to the 
low frequency one). 

Changing continuously the feedback phase by ±27r 
transforms a CS solution into a different one, associated 
to an adjacent longitudinal external cavity mode. For ex- 
ample, increasing from to 27r transforms Z\ into Z3. 
We cannot speak here of multistability because most of 
these resting solitonic states are unstable for the param- 
eters we have considered. We will show below, however, 
that these unstable states can be stabilized by the effect 
of the carrier dynamics in model ([1]), leading to bistabil- 
ity and apparently to multistability. 
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FIG. 8: Dynamics of the system ([2]) for an initial condition 
just above the point Z2 of Fig. (^) Time-space plot, 

where the vertical direction is x and the horizontal direction 
is time, on the same scale as in panel (b). White corresponds 
to large amplitude, black to zero, (b) Time evolution of the 
CS maximum amplitude 
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FIG. 9: Dynamics of the system, plotted as in Fig. [8] for an 
initial condition just above the point Za of Fig. [7{a). 

Fig. [7] also shows the branch of self-propelled CSs. 
Open (filled) circles indicate unstable (stable) solutions. 
For simplicity we have not plotted here the equivalent 
solutions obtained for a 27r shift in the feedback phase. 
As in the case of resting solitons, there can be multi- 
ple moving states and, moreover, they can be bistable 
because the interval of feedback phases where they are 
stable is slightly broader than 27r. 

Finally Figs. [5] and [SI show two examples of the dy- 
namics of the system starting from an initial condition 
just above the lower unstable CS branch in Fig.[7l^a). 

In the first case (Fig. [5]) we start the simulation with an 
initial condition just above point Z2 in Fig.[7]Ja). We ob- 
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serve first a slow evolution away from the lower unstable 
CS to form a transient approximation to the unstable 
state Zi. We note that Zi differs from initial state in 
both amplitude and frequency. Then the amplitude and 
frequency switch to those of the stable self-propelled CS 
corresponding to the filled (red) circle intersecting the 
vertical dashed line (/) = in Fig. [71(a) . Such evolution 
is typical for feedback phase values for which the only 
stable CS state is the self-propelled CS, and for initial 
conditions above the lower branch solitonic states. 

There is a second scenario, typical for values of the 
feedback phase for which both the resting and moving 
upper-branch CSs are unstable. An example is shown 
in Fig. [9l Starting just above the point the system 
approaches Z^^ and then starts oscillating with a large 
amplitude before switching off completely. This long 
transient excursion, for a perturbation just above the CS 
lower branch is an indication of excitable behavior, since 
perturbations just below this threshold decay directly to 
the zero solution. The period of the oscillations is close 
to T, as might be expected, which indicates that more 
than one longitudinal mode is excited. Initially, it seems 
that there is a beating of the two CS modes Z^ and Zi, 
but in the later stages the sharpness of the spikes (both 
dark and, latterly, bright) indicates transient locking of 
at least three CS modes. 



V. INFLUENCE OF CARRIER DYNAMICS 

In this section we address the influence of the carrier 
lifetime on the dynamics of the system, which is typi- 
cally very important in semiconductor laser media. For 
this, we use Eqs. (P) instead of the simplified version ([2]), 
studied in previous sections. The resting solitons are ac- 
tually solutions of both models ^ and ([T|) , because the 
adiabatic elimination of the carriers does not influence 
these steady states. As a prototypical example to study 
the effects of the carrier dynamics on the CS, we will an- 
alyze how the stability of state Zi of Fig. [7] changes with 
the carrier relaxation time Ti . 

For small enough values Ti < Tf 0.085 ± 0.001, the 
behavior of Eqs. ([T|), starting from the solitonic state Zi, 
differs from the behavior of Eqs. ^ only quantitatively, 
i.e. we observe in both cases the transition from a resting 
to a self-propelled CS. Fig. [10] shows this dynamics for 
Ti = 0.05, which can be compared with the middle and 
final part of Fig. [5] where a soliton started from Z2 first 
evolves towards Zi and then starts to move. 

For Ti > Tf , the drifting instability is transformed 
into a swinging instability, where the CS starts to per- 
form spatial oscillations around its initial position (see 
Fig. [TTk). There are several physical mechanisms that 
may contribute to this behavior. One is delay, that in 
combination with the motion of the structure gives a pos- 
itive feedback at the rear part of the CS, slowing down 
the motion, and possibly causing the oscillations. An- 
other one is the slow relaxation of the carriers, which 



(a) 




t 



FIG. 10: Dynamics of system ([T]), plotted as in Fig.[8l starting 
from an initial condition corresponding to the point Z\ of 
Fig. [7] Here Ti = 0.05. 



(a) 




t 



FIG. 11: The same as in Fig. [TO]for Ti = 0.5. 



also opposes the unidirectional movement of the CS. The 
oscillations in space arc also accompanied by intensity 
oscillations of the maximum fFig. [TTb i. The amplitude 
of the oscillations grows initially, then evolves in a com- 
plex manner before the CS finally switches off. The full 
quasi-period of the oscillations is again close to r, indi- 
cating multi-mode behavior, though there are two "max- 
imum maxima" per period, and hence two oscillations 
in Fig. Illf b) per feedback time. Swinging, followed by 
switch off, is observed for an interval of carrier relaxation 
times between and T'^ ^ 0.791 ± 0.001. 

For carrier relaxation times larger than the state Z5 
becomes stable, and the swinging instability of the state 
Zi leads to the birth of a pair of solitonic states Z5, as 
shown in Fig. [Hj Finally, for Ti > Tf = 1.885 ± 0.005 
the state Zi is stable. Moreover, state Z5 remains stable, 
so therefore in this region there is bistability between two 
resting CSs with different amplitude, which is not found 
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FIG. 12: (Color online). Dynamics of the system with the 
same parameters as in figures [10] and 1111 except Ti = 1. 
(a) Space-time plot of the field amplitude for time interval 
[750;1000]. (b) Evolution of the CS maximum amplitude; (c) 
Location of the CS majdmum, the split signals the birth of 
two states; (d) Evolution of the instantaneous CS frequency 
u) = Arg{Emax{t) / Emax{t — At))/At, all for time interval 
[0;3500]. 



in dH). 

For the chosen value of the detuning lJs, and typical 
values of Ti, self-propelled self-localized states disappear. 
However increasing the detuning, which reinforces the 
drifting instabilities as discussed in Sec. Illli we observe 
stable self-propelled ID CSs for Eqs. ^ for Ti = 1, 
ujs = 200, fjL = 0.8, with the other parameters as in Fig. [51 
We note that this value of Ti is typical [H, [l^l for carrier 
relaxation in VCSELs. 

Summarizing this Section, carrier dynamics plays a re- 
ally important role in the evolution of the system. For 
relatively long carrier lifetimes it allows bistability of rest- 
ing CSs in parameter regimes where it is not present in 
0. In an intermediate range of lifetimes it can lead to 
a new "swinging" CS instability, which in turn may lead 
to "CS fission" . 



VI. FINAL REMARKS 



We have studied in detail the formation of self- 
propelled localized solutions in a Cavity Soliton Laser 
composed of a VCSEL subject to filtered external optical 
feedback. These states are potentially useful for applica- 
tions such as all-optical delay lines p. Our results have 
been obtained for systems with one spatial dimension, 
but evidence of a qualitative agreement with the more 
realistic 2D case is given. 

The self-propelled soliton branch has been constructed 
by solving with a Newton method the stationary equa- 
tions in the reference frame moving with the soliton. The 
branch originates from the modulational instability bi- 
furcation point of the trivial solution where tilted waves 
appear with the critical wavenumber. Apparently, ampli- 
tude modulation of the tilted wave leads to the existence 
of self-propelled CSs in such systems. 

We have also shown that the characteristics of CSs can 
be controlled via the feedback phase, and the existence 
of multiple localized solutions associated to different ex- 
ternal cavity modes. 

We have analyzed the influence of the carrier relax- 
ation time on the properties of the solitonic states. It 
has been shown that increasing Ti leads to the stabiliza- 
tion of resting solitonic states. This opens the possibility 
of observing multistability of CSs. For intermediate val- 
ues of the carrier relaxation time we have observed a new 
kind of instability leading to the oscillation of the posi- 
tion of the CS and, for larger Ti, to the formation of two 
identical solitons at different positions. 

Finally, we should mention, that the models used here 
still represent only a rather rough approximation to ex- 
periment. The complexity of real devices is much higher, 
involving different spatial mechanisms such as carrier dif- 
fusion [2l| . gain and loss dispersion [2^ . thermal fre- 
quency shi ft I23| and multiple round-trips in the exter- 
nal cavity [24 l25l. [26j . Simple models such as ours are, 
nonetheless, useful for understanding fundamental fea- 
tures of basic phenomena such as self-propelled solitons, 
and for the prediction of some new phenomena, like the 
swinging instabilities. 
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